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ABSTRACT 

We present a linear analysis of the vertical structure and growth of the magnetoro- 
tational instability in stratified, weakly ionised accretion discs, such as protostellar 
and quiescent dwarf novae systems. The method includes the effects of the magnetic 
coupling, the conductivity regime of the fluid and the strength of the magnetic field, 
which is initially vertical. The conductivity is treated as a tensor and assumed constant 
with height. 

We obtained solutions for the structure and growth rate of global unstable modes 
for different conductivity regimes, strengths of the initial magnetic field and cou- 
pling between ionised and neutral components of the fluid. The envelopes of short- 
wavelength perturbations are determined by the action of competing local growth 
rates at different heights, driven by the vertical stratification of the disc. Ambipolar 
diffusion perturbations peak consistently higher above the midplane than modes in- 
cluding Hall conductivity. For weak coupling, perturbations including the Hall effect 
grow faster and act over a more extended cross-section of the disc than those obtained 
using the ambipolar diffusion approximation. 

Finally, we derived an approximate criterion for when Hall diffusion determines 
the growth of the magnetorotational instability. This is satisfied over a wide range of 
radii in protostellar discs, reducing the extent of the magnetic 'dead zone'. Even if 
the magnetic coupling is weak, significant accretion may occur close to the midplane, 
rather than in the surface regions of weakly-ionised discs. 

Key words: accretion, accretion discs - instabilities - magnetohydrodynamics - 
stars: formation. 



1 INTRODUCTION 

The collapse of protostellar cores lead to the development of 
a central mass or protostar, surrounded by a disc of material, 
which is accreted towards the center. During this process 
angular momentum is transferred to a small percentage of 
disc material at large radii, enabling the collapse of most of 
the disc towards the central star (e.g. Weintraub, Sandell 
& Duncan 1989; Adams, Emerson & Fuller 1990; Beckwith 
et al. 1990). The evolution of this 'disc accretion' phase is 
dependent upon the rate of angular momentum transport in 
the disc (e.g. Adams & Lin 1993). 

A variety of mechanisms have been invoked to explain 
this transport. As the molecular viscosity of accretion discs 
is too low to explain observed accretion rates (Pringle 1981), 
some form of turbulent viscosity must be present. The origin 
and characteristics of this turbulence remains an important 
problem in star formation theories. Convective turbulence 
has been considered as an option (Lin & Papaloizou 1980), 
but further studies suggest that this mechanism may trans- 



port angular momentum towards the central star instead 
of away from it (Cabot & Pollack 1992; Ryu & Goodman 
1992; Stone & Balbus 1996). The gravitational field of a 
companion star may trigger hydrodynamic waves which can 
transport angular momentum (Vishniac & Diamond 1989; 
Rozyczka & Spruit 1993), but as a significant fraction of 
stars do not belong to binary systems, this mechanism is not 
general enough to explain accretion processes in all stars. 

Balbus and Hawley have pointed out that the nature 
of this anomalous viscosity can be hydromagnetic (Bal- 
bus & Hawley 1991; Hawley & Balbus 1991; Stone et al. 
1996). This 'magnetorotational' instability (MRI) had been 
described initially by Velikov (1959) and Chandrasekhar 
(1961) through their analysis of magnetised Couette flows. 
It drives turbulent motions which transport angular momen- 
tum radially outwards, as fluid elements exchange angular 
momentum non-locally by means of the distortion of the 
magnetic fleld lines that connect them. 

Under ideal MHD conditions, MRI perturbations grow 
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in discs that are differentially rotating, with the angular 
velocity increasing outwards. Axisymmetric modes need a 
magnetic field with a weak, poloidal component. In this con- 
text, 'weak' means that the magnetic energy density of the 
field is less than the thermal energy density. These pertur- 
bations have a characteristic length scale A ~ v^/O, where 
va is the Alfven speed and f2 is the keplerian angular fre- 
quency in the disc, and a maximum growth rate ~ 9^2/2, 
with q = 1.5 for keplerian discs (Balbus & Hawley 1992a). 
This growth rate docs not depend on the strength or direc- 
tion of the magnetic field as long as a poloidal component is 
present. Non-axisymmetric perturbations are most unstable 
under the influence of a poloidal field, but also grow at a 
reduced rate if the field is purely toroidal (Balbus & Hawley 
1992b). These perturbations are of interest for the analysis 
of field amplification mechanisms, as dynamo amplification 
can not occur through axisymmetric perturbations (Moffatt 
1978). With no strong dissipation processes, no other condi- 
tions are required. Because of its robustness and the general 
conditions under which it develops, the rnagnctorotational 
instability is a promising source of turbulent viscosity in ac- 
cretion discs. 

Ideal MHD conditions are a good approximation to 
model astrophysical systems where the ionisation fraction 
of the gas is high enough to ensure neutral and ionised com- 
ponents of the fluid are well coupled. Active dwarf novae 
(except possibly in the outer regions) and black hole accre- 
tion discs arc examples of such systems. However, in dense, 
cool environments such as those of protostellar discs, it is 
doubtful that magnetic coupling is significant over the en- 
tire radial and vertical dimensions of the discs (Gammie 
1996; Wardle 1997). Similar conditions are thought to apply 
in quiescent and the outer regions of hot-state dwarf novae 
discs (Gammie & Menou 1998; Menou 2000; Stone et al. 
2000). In these cases, low conductivity significantly affects 
the growth and structure of MRI perturbations. Different 
approximations have been adopted to account for the de- 
parture from ideal MHD in low conductivity astrophysical 
discs (see section 2.1). 

Most models of the MRI in low conductivity discs have 
used the ambipolar diffusion (Blaes & Balbus 1994; MacLow 
et al. 1995; Hawley & Stone 1998) or resistive (Jin 1996; 
Balbus & Hawley 1998; Sano, Inutsuka & Miyama 1998) 
limits. Recently, it has been recognised the importance of 
the Hall conductivity terms in addition to resistivity for the 
analysis of low conductivity discs (Wardle 1999, hereafter 
W99, Balbus & Terquem 2001 and Sano & Stone 2002a,b; 
2003). 

The huge variation of fluid variables over the vertical 
and radial extension of astrophysical discs is a further com- 
plication. Vertical stratification is particularly relevant, as 
these objects are generally thin and changes in the plane 
of the disc are much more gradual than those in the direc- 
tion perpendicular to it. Previous models of the MRI have 
not included density stratification and Hall conductivity si- 
multaneously. It is expected that solutions will be strongly 
modified when both factors are present. This motivates the 
present study. 

This paper examines the structure and linear growth 
of the magnetorotational instability in vertically stratified, 
non self-gravitating accretion discs. We assume the disc is 
isothermal and geometrically thin, so variations in the fluid 



variables in the radial direction can be ignored. The initial 
magnetic field is vertical and the analysis is restricted to 
perturbations with wavevector perpendicular to the plane 
of the disc (fc = fcz). These are the most unstable perturba- 
tions with the adopted field geometry, as magnetic pressure 
strongly supresses displacements with kr ^ (Balbus & 
Hawley 1991, Sano & Miyama 1999). The conductivity of 
the gas is treated as a tensor and assumed constant with 
height in this initial study, although the formulation is also 
valid for a z-depcndcnt conductivity. This makes the present 
method a powerful tool for the analysis of more realistic discs 
(see section 6 - Discussion). 

Section 2 presents the governing equations for a weakly 
ionised, magnetised disc in ncar-keplcrian motion around the 
central star and details the adopted disc model. Section 3 
summarises the linearisation of the equations and presents 
the final linear system in dimensionless form. It also de- 
scribes the three parameters that control the dynamics of 
the fluid. Section 4 discusses the boundary conditions used 
to integrate the equations from the midplane to the surface 
of the disc and the integration method. Section 5 presents 
the test cases used to characterise the conductivity regimes 
relevant for this work and compares our results with a previ- 
ous local analysis. It also details key findings on the depen- 
dency of the structure and growth rate of the perturbations 
with the conductivity regime, the strength of the magnetic 
field and its coupling with the neutral gas. These results 
arc discussed in section 6. By way of example, this section 
also calculates the structure and growth rate of the MRI 
under fluid conditions where different conductivity regimes 
are dominant at different heights above the midplane. This 
reflects (qualitatively) the conditions expected to be found 
in real discs. Finally, the methodology and key findings of 
this paper are summarised in section 7. 



2 FORMULATION 

2.1 Governing Equations 

The equations of non-ideal MHD are written about a lo- 
cal keplerian frame corotating with the disc at the angular 
frequency f2 associated with the particular radius of inter- 
est. Consequently, the velocity of the fluid can be expressed 
as a departure from exact keplerian motion v = V — vk, 
where V is the velocity in the standard laboratory coordi- 
nate system (r, ^, z) anchored at the central mass M, and 
Vk = \/ GM/r4> is the keplerian velocity at the radius r. 
Similarly, if d/dt is the time derivative in the local keple- 
rian frame, then the time derivative in the laboratory frame 
can be expressed as d/dt + ild/dcp. We also assume that 
the fluid is weakly ionised, meaning that the abundances 
of charged species are so low that their inertia and thermal 
pressure, as well as the effect of ionisation and recombination 
processes in the neutral gas, are negligible. These assump- 
tions effectively restrict the range of frequencies that can be 
studied with this formulation to be smaller than the collision 
frequency of any of the charged species with the neutrals. 
Accordingly, separate equations of motion for the charged 
species are not required and their effect on the neutrals is 
contained in a conductivity tensor (see section 2.2). 

The governing equations are the continuity equation, 
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^+V.(p.)=0, 
the equation of motion, 
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and the induction equation, 
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In the equation of motion (2), $ is the gravitational 
potential due to the central gas, given by 



$ = -- 



GM 



(r2 + 02)2 



(4) 



and v\lr is the centripetal term generated by exact kep- 
lerian motion. At the disc midplane this term balances the 
radial component of the gravitational potential. The terms 
2Qv^f and ^Qvr4> are the coriolis terms associated with the 
use of a local keplerian frame, Cs = \pPjp is the isothermal 
sound speed, Q. = vk It is the keplerian frequency and c is 
the speed of light. Other symbols have their usual meanings. 

In the induction equation (3), the term cV X iJ'contains 
the effects of non- ideal MHD. E' is the electric field in the 
frame comoving with the neutrals and the term |QBr<^ ac- 
counts for the generation of toroidal field from the radial 
component due to the differential rotation of the disc. 

Additionally, the magnetic field must satisfy the con- 
straint: 



V-B =0, 

and the current density must satisfy Ampere's law, 



J= —V X B 

4-w 



and Ohm's law, 
J = a E' . 



(5) 
(6) 
(7) 



Note that the conductivity, which depends on the abun- 
dance and drifts of the charged species through the neutral 
gas is treated as a tensor cr, as detailed in the following 
section. This formulation is compared to the drift velocity 
approach in section 2.3. 



2.2 The conductivity tensor <t 

The electric conductivity is a tensor whenever the gyrofre- 
quency of the charged carriers is larger than the frequency 
of momentum exchange by collisions with the neutrals, or 
\(3\ > 1 (Cowling 1957, Norman & Hayvaerts 1985, Nakano 
& Umebayashi 1986) . On the contrary, when collisions with 
the neutrals are dominant the conductivity is a scalar, the 
ordinary ohmic resistivity. To obtain expressions for the 
components of this tensor, we begin by writing down the 
equations of motion of the ionised species. As inertia and 
thermal pressure are neglected, the motion of the charged 
particles is given by the balance of the Lorentz force and the 
drag force from collisions with the neutrals, 



Zie(E' + 



X B 



) - ^jmjVj = , 



(8) 



where each charged species j is characterised by its number 
density rij, particle mass rrij, charge ZjC and drift velocity 
Vj . In the above equation, 

< av >j 



rrij + m 



(9) 



where m is the mean mass of the neutral particles and 

< av >j is the rate coefficient of momentum exchange by 
collisions with the neutrals. We will also make use of the 
Hall parameter, 

Z,eB 1 



(10) 



given by the ratio of the gyrofroquency and the collision fre- 
quency of charged species j with the neutrals. It represents 
the relative importance of the Lorentz and drag terms in 
equation (8). 

Following the treatment of Wardlo & Ng (1999) and 
W99 we use the following expression for Ohm's law. 



J = (T E' = a^En + aiB X E'j_ + a2E'j_ , 



(11) 



obtained by inverting (8) to express Vj as a function of 
E' and B and then using J = e'}2j ^ijZjVj together with 
the charge neutrality assumption Y2j ^jZj = 0. In equation 
(11), and E'j_ are the components of the electric field 
E' , parallel and perpendicular to the magnetic field, respec- 
tively. The components of the conductivity tensor <t are the 
conductivity parallel to the magnetic field, 
ec ■ 



(12) 



the Hall conductivity, 
ec ^ rijZ 



(71 



31' 



and the Pedersen conductivity. 



CT2 



B^ 1 -H 



(13) 



(14) 



The relative values of the components of the conductiv- 
ity tensor differentiate three conductivity regimes: 

(i) The ambipolar diffusion regime occurs when cr|| ^ 
o"2 3> |<Ti|, or ^ 1 for most charged species. This im- 
plies that most charged particles are strongly tied to the 
magnetic field by electromagnetic stresses. This regime is 
dominant at relatively low densities, where the magnetic 
field is frozen into the ionised component of the fluid and 
drifts with it through the neutrals. The linear behaviour of 
the MRI in this regime has been analysed by Blaes & Balbus 
(1994) and the non-linear growth by MacLow et al. (1995) 
and Hawloy & Stone (1998). 

(ii) The resistive (Ohmic) regime is obtained when most 
charged species are linked to the neutrals via collisions. This 
occurs when a\\ w (72 S> \ai\, implying <C 1. This regime 
is predominant closest to the midplane, where the high den- 
sity makes the collision frequency of the charged particles 
with the neutrals high enough to prevent the former from 
drifting. This case has been studied under a linear approxi- 
mation by Jin (1996), Balbus & Hawley (1998), Papaloizou 
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& Terquem (1997), Sano & Miyama (1999) and Sano et. 
al. (2000). The non- linear regime has been studied by Sano 
et al. (1998); Fleming, Stone & Hawley (2000) and Stone 
& Fleming, (2003). This last work includes a a-dependent 

resistivity. 

(iii) Finally, the Hall regime occurs when charged parti- 
cles of one sign axe tied to the magnetic field while those of 
the other sign follow the neutrals. In this case \ai \ ~ (72 < cr\\ 
and \f3\ fs 1. It is important at intermediate densities, be- 
tween the ones associated with the ambipolar and Ohrnic di- 
fussion regimes. Recent studies have explored the MRI with 
Hall effects in the linear (W99, Balbus & Terquem 2001) 
and non- linear regimes (Sano & Stone 2002a,b; 2003). 

2.3 Comparison with the multifluid approach 

Another commonly used form of the induction equation is 
obtained by assuming that ions and electrons are the main 
charge carriers and drift through the neutrals (e.g Balbus 
& Terquem 2001, Sano & Stone 2002). Using this approach, 
the induction equation is 



^=Vx(.. 
= V X 



X B) 



vx B- 



J X B ^{J X B)x B 



CJiPPi 



(15) 

where Ve is the electron drift speed, r] = c^/A-na is the re- 
sistivity and subscripts e and i refer to electrons and ions, 
respectively. The four terms in the right hand side of equa- 
tion (15) are, from left to right, the inductive, resistive. Hall 
and ambipolar diffusion terms. We now express equation (3) 
in terms of J, B and the components of the conductivity 
tensor in order to show that, in the appropriate limit, it 
corresponds to equation (15), as expected. Wo begin by in- 
verting equation (11) to find an expression for E' , 



, J ai J X B 

^ = 1 2 D 



02 



1 



{J X B)xB 



(16) 



where crx = -^Z af + is the total conductivity perpendic- 
ular to the magnetic field. Assuming that the only charged 
species are ions and electrons with Hall parameters /3i and 
0e (< 0). respectively, and that charge neutrality is satis- 
fied (rii = rie), we obtain the following expressions for the 
components of the conductivity tensor: 
eerie 



5 -(A -A), 

cene il3i + I3e){l3e - Pi) 
B (i + /3|)(i + /32) 

cen, {I ~ (3if3,){Pi - f3e] 
B (l + /3|)(l + /32) 

From (18) and (19), we find: 

(A - /3e) 



CT2 



and 



(T± = 



B [(i + /3|)(i + /32)]i/2- 
Substituting these expressions into (16) gives 



(17) 
(18) 
(19) 

(20) 



E' = 



J ^ /3e + A J X B 

(T Pe— Pi CeUe Pe — A 



Pe {J X B)x B 



(21) 

where a is the electrical conductivity duo to electrons. Fi- 
nally, simplifying (21) by using |/3e| ^ Pi (as is the case) 
and substituting the resulting E' in equation (3) (without 
the coriolis term), yields the standard result, shown in (15). 

Each of the last three terms in the right hand side of 
equation (15) dominate when the fiuid is in a particular 
conductivity regime (section 2.2). These limits can also be 
recovered with the appropriate assumptions, through equa- 
tions (16) and (21). To get the resistive regime, for example, 
we substitute ai = and ay — a-i into equation (16). In this 
limit the conductivity is a scalar (the resistive approxima- 
tion), so E' = J I a and the induction equation reduces to 
the familiar form. 



dB 

dt 



V X 



V xB ■ 



4na 



V X B 



(22) 



On the other hand, to model the Hall limit, we regard 

the charged species to be either 'ions', which are strongly 
tied to the neutrals through collisions (Pi <C 1), or 'elec- 
trons', for which the only important forces are electromag- 
netic stresses {\pe\ ^ 1) (see also discussion in W99). In this 
limit the Hall conductivity is 



B 



while the Pedersen conductivity is 



(T2 



0-1 



A-^l«ar. 



(23) 



(24) 



As the ions are effectively locked with the neutrals, the 
current density will be given by the drift of the electrons 
through the neutral gas. Collisions are unimportant in their 
equation of motion, so they drift perpendicular to the plane 
of the electric and magnetic fields, in order to annul the 
Lorentz force acting upon them. In this limit {J X B) X B = 
—J and 



E' = 



J X B 



eerie, 



(25) 



consistent with the Hall term in equation (15). 

Filially, tiie ambipolar diffusion approximation is recov- 
ered by assuming a\\ 3> (72, ai = and |/3e| ^ /3i S> 1. In 
this limit, the Pedersen conductivity is given by 



<72 



ce Ui 

Bji' ^""^ 



E' = - 



(J xB)xB 

c'^-jiPPi 



(26) 



(27) 



which originates the ambipolar diffusion term of (15). 

Although the multifluid drift and conductivity tensor 
formulations are ultimately equivalent, which one is more 
convenient depends on the problem at hand. In particular, 
the presence of dust grains tends to make the treatment 
of different species especially complex. In protostellar discs, 
dust grains can be the more abundant charged species over 
extended regions. For example, assuming 0.1 p,m grains, neg- 
atively charged grains dominate whenever uh ^ 10^^ cm~^, 
while positive charged grains are the most abundant ions 
for UH > 10" cm"^ (Wardle & Ng 1999). Having separate 
equations of motion for different charged species would gen- 
erally involve dividing the grain size distribution of interest 



Magnetorotational instability in stratified, weakly ionised accretion discs 5 



into an appropriate number of discrete intervals and explic- 
itly treating each one. Unless the number of such intervals 
is small it is easy to see that this method can become very 
cumbersome. In these circumstances, incorporating the con- 
tribution of each charged species into a conductivity tensor 
can be a valuable approach. 

2.4 Disc Model 

Our model incorporates the vertical structure of the disc, 
but neglects fluid variations in the radial direction. This 
is appropriate as astrophysical accretion discs are generally 
thin and changes in the radial direction occur in a much 
bigger length scale than those in the vertical direction. In- 
cluding the vertical structure means that perturbations of 
spatial dimensions comparable to the scale height of the disc, 
which are associated with a strong magnetic field {va ~ c^), 
or low conductivity, can be explored. 

The balance between the vertical component of the cen- 
tral gravitational force and the pressure gradient within the 
disc determines its equilibrium structure. The vertical den- 
sity distribution in hydrostatic equilibrium is given by 

In the above equation, Po{r) is the gas density at the mid- 
plane and H{r) = Cs/fi is the scaleheight of the disc. 

A self-consistent treatment of this problem, would in- 
volve adopting a particular dependency of po and H with 
r using a suitable model, such as the minimum solar neb- 
ula (Hayashi, Nakazawa & Nakagawa 1985) and calculating 
p{r,z) by means of (28). This density, together with the 
adopted strength of the magnetic field B and the values of 
the conductivity tensor (t as a function of height would be 
used to evaluate the parameters that govern the fluid evo- 
lution (see section 3.3) and solve the fluid equations. 

The realistic evaluation of the conductivity tensor is a 
complex undertaking, as it depends critically on the abun- 
dances of charged species (ions, electrons and charged dust 
grains) which, in turn, are a function of the ionisation bal- 
ance in the disc. This balance is given by the equilibrium 
between ionisation processes by cosmic rays, radioactive el- 
ements and X-rays from the central star (e.g. Hayashi 1981; 
Glassgold, Najita & Igea 1997, Igea & Glassgold 1999; Fro- 
mang, Terquem & Balbus 2002) and recombinations taking 
place both in the gas phase and in grain surfaces (e.g. Nishi, 
Nakano & Umebayashi 1991). In low conductivity discs, the 
level of ionisation is insufficient to produce good coupling 
between the magnetic field and the neutral component of 
the fluid over their entire vertical structure. In protostcUar 
discs, for example, it is expected that in the region outside 
~ 0.1 AU from the central star, the coupling will be sig- 
nificant only in the surface layers, where X-rays and cosmic 
rays can penetrate and ionise the gas (Gammie 1996; Wardle 
1997). In these environments, the ^-dependent atenuation of 
the ionisation rate typical of interstellar space 10~^'^H~^ 
has to be taken into account. 

On the other hand, the contribution of dust grains to re- 
combination processes is particularly complex, because they 
generally have a distribution of sizes and corresponding colli- 
sion cross sections (Mathis, Rumpl & Nordsieck 1977; Ume- 
bayashi & Nakano 1990). Moreover, the dynamics of grain 



particles depends on the activity of the disc. In quiescent 
environments, they tend to settle towards the midplane and 
begin to agglomerate into bigger structures that could even- 
tually become planets (e.g. Weidenschilling & Cuzzi 1993). 
This removes grains from regions at relatively high z and 
causes the ionisation level to increase by reducing recombi- 
nation processes taking place in their surfaces. Simulations 
of the evolution of dust grain distributions in dense cores 
show that they grow icy mantles and coagulate efficiently 
(Osscnkopf 1993). By this process, the smallest particles 
grow quickly while the upper grain size limit is only slightly 
changed (Ossenkopf 1993). As a result, the surface area of 
dust grains can be significantly modified by grain evolution, 
which ultimately affects the ionisation balance in the disc. 

In the present work, a simpler treatment has been 
adopted. In order to study the MRI under different con- 
ductivity regimes, the values of the components of the con- 
ductivity tensor have been selected so that they satisfy the 
conditions outlined in section 2.2 for each regime. For sim- 
plicity, those values are assumed to bo constant, although 
the formulation allows them to be a function of height. We 
found that the parameters that control the evolution of the 
fiuid are the ratio of the components of the conductivity ten- 
sor perpendicular to the magnetic field {ai/a2), the strength 
of the magnetic field and its degree of coupling with the 
neutral component of the fluid (see section 3.3). The mid- 
plane values of these parameters have been selected in order 
to simulate the fluid conditions we wore interested in mod- 
elling. This approach will be discussed in section 5.1 where 
the chosen test cases axe detailed. 

With this approach we are able to study the depen- 
dency of the growth rate and structure (characterised by 
the height of maximum amplitude and the wavermmbor of 
unstable modes) of the instability with the parameters of 
the fluid in a stratified disc. This is relevant, as the re- 
gion where linear perturbations peak is also cxpcc:tcd to be 
the region where non-linear perturbations grow fastest, un- 
til turbulence finally sets in and causes all wavenumbers to 
interact so that, eventually, the longest wavelengths carry 
the greatest angular momentum transport (e.g. Hawley & 
Balbus 1995). 



3 LINEARISATION 

We linearised the system of equations (1) - (3), (6) and (11) 
about an initial steady state where J = v = E' = Q and 
B = Bz. In the initial state both E' and J vanish, so the 
changes in the conductivity tensor due to the perturbations 
do not appear in the linearised equations. As a result, it is 
not necessary to explore how the perturbations affect the 
conductivity and only the values in the initial steady state 
are required. 

3.1 Linearised Equations 

We assume the wavevector of the perturbations is perpen- 
dicular to the plane of the disc {k — kz). Perturbations with 
a vertical wavermmbor, initiated from a vertically aligned 
equilibrium magnetic field, exhibit the fastest growth rate 
for a given set of parameters because magnetic pressure 
strongly suppresses displacements with kr =^ (Balbus & 
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Hawley 1991; Sano & Miyama 1999). Taking perturbations 
of the form q = Qq + 6q{z)e'''^^ about the initial state, lin- 
earising and neglecting terms of order H/r or smaller, we 
find that the equations decouple into two subsystems. One 
of them corresponds to sound waves propagating in the vor- 
tical direction and the other describes MHD perturbations 
in the plane of the disc, with vanishing z component. 

With these simplifying assumptions, the final linear sys- 
tem of equations that describes the MHD perturbations 
within the disc is. 



iujpdvr — 2pQ5v^ dJ,j, = , 

c 



0, 



iujSBr — c 



iuiSBA, + c 



SJr 



dSE4 
dz 

dSEr 

dz 
c dSB, 
4n dz 
c ddBr 
A-K dz 



= 0, 



+ ^Q.5Br=Q, 



(29) 
(30) 
(31) 
(32) 
(33) 
(34) 

(35) 
(36) 



where SE^ and SEr arc the perturbations of the electric 
field in the laboratory frame, given by 



6E^ = 5E'^ H 5vr , and 



6Er = SE' - —5vi 



(37) 
(38) 



We note that a\\, the component of the conductivity 
tensor parallel to the magnetic field, does not appear in the 
final linearised equations. Because the ambipolar diffusion 
and resistive conductivity regimes are differentiated by the 
value of cr|| (section 2.2), this means that during the linear 
stage of the MRI, under the adopted approximations, these 
regimes are identical. 

3.2 Equations in dimensionless form 

Equations (29) to (38) can be expressed in dimensionless 
form, normalising the variables as follows: 



z 
H 



P = 



p{r,z) 
Po(r) 



SB* = 



5B 

Bo 



6v* = 



Sv 



SE* = 



cSE 
CsBo 



SE'* 



cSE' 
CsBo 



5 J* = 



cSJ 



cr± = 



cr±o 



CsBo(T±g a±Q 

Subscript 'o' is used to denote variables at the mid- 
plane of the disc. Effecting these changes and dropping the 
asterisks to keep the notation simple, we finally express the 
dimensionless system of equations in matrix form as: 





( Br \ 


1 








CiAi 


CIA2 


d 












-C1A2 


CiA-, 


dz 






3 
2 


—l> 












V 


V 












/ Br \ 

\E^ ) 



(39) 



SJ = a±C2 



Sv = 



Xo 



A2 -A3 



Al A2 



-2 



SE 



p l + 



SJ 



SE' = 
where 



(72 0"l 



— fTi —(72 



SJ 



Ci = XoCr± ( ^ 1 C2; 



VA 



C2 



^1 



A2 = 



1 + 



(Tl 

02 
a± 

(Tl 

(T± 



+ 2 



+ 1^ 



P 

Xo<T± 



1^+2.^ + ^ 
1 + v' \2 a± a± p 



1 



p 


1 + 1/2 


Xo(T± 


1 


p 


l + Z/2 


Xo<T± 


1 



and 



2 p H- j/2 
In the above expressions. 

Bo 



VA 



^JA-KPo 

is the Alfven speed at the midplane of the disc, and 



Xo = 



Oco _ 1 B„a±o 



(40) 

(41) 

(42) 

(43) 
(44) 

,(45) 
(46) 
(47) 
(48) 

(49) 
(50) 



is a parameter that characterises the midplane coupling be- 
tween the magnetic field and the disc (see section 3.3). To 
understand the information contained in this parameter it 
is useful to recall that the effect of finite conductivity is 
different for perturbations of different wavelengths. Finite 
conductivity is important when the term cV X E' in the 
induction equation (3) is of order V X (t; X B). Adoptinga 
length scale L ~ 1/k ~ va/u> and a --^ ax, it is found that 
these two terms are comparable when kvA ~ (c^/47r)fe^/crx, 
in other words, non-ideal effects will strongly modify wave- 
modes at, or above, the critical frequency. 



(51) 



Magnetorotational instability in stratified, weakly ionised accretion discs 7 



It can be shown (W99) that in the limit |/3jj oo, ujc 
reduces to Y2j^jPjy ^^e collision frequency of the neutrals 
with any of the charged species. Generally, uJc is smaller 
than this value and much smaller than -yjp, the collision 
frequency of charged species j with the neutrals. In dense 
clouds, LOc ~ 7gPg, which is the smallest 7jPj, the col- 
lision frequency of neutrals with grains. As pa ~ O.Olp, 
the treatment of this paper, restricted to w < 7jp by ne- 
glecting the inertia of the charged species, remains valid for 
Lo ~ lOOcJc. For perturbations with lower frequencies (longer 
wavelength) than ujc, ideal MHD (the flux- freezing approx- 
imation) is valid. 



3.3 Parameters 

As these equations reveal, three important parameters con- 
trol the evolution of the fluid: 

(i) va/cs, the ratio of the Alfvon speed to the isothermal 
sound speed of the gas at the midplanc. It is a measure of the 
strength of the magnetic field. In ideal MHD unstable modes 
grow when the magnetic field is subthermal {va/cs < 1). 
When va ~ Cs the minimum wavelength of the instability 
is of the order of the scaleheight of the disc and the growth 
rate decreases rapidly. 

(ii) Xo, a parameter that characterises the strength of the 
coupling between the magnetic field and the disc at the mid- 
plane (see equation 50) . It is given by the ratio of the critical 
frequency above which flux-freezing conditions break down 
and the dynamical frequency of the disc at the midplane. If 
Xo = ojco/^ < 1 the disc is poorly coupled to the disc at the 
frequencies of interest for dynamical analysis. As the growth 
rate of the most unstable modes are of order Q, in ideal MHD 
conditions, these are also the interesting frequencies for the 
study of this instability. 

(iii) (Ti /a2 , the ratio of the conductivity terms perpendic- 
ular to the magnetic field. It is an indication of the conduc- 
tivity regime of the fluid, as discussed in section 2.2. 

Note that the density of the disc decreases with z, so 
the local values of x and va/cs increase with height. The 
parameters of the model are defined as the corresponding 
values at the midplane. 

It is common practice to characterise the magnetic cou- 
pling of a weakly ionised fluid by its electron density rie . Be- 
fore finishing this section, we discuss how x relates to this 
fluid parameter. We begin by writing the magnetic Reynolds 
number as (e.g. Balbus & Terquem 2001), 



ReM = 



vaH 
V 



a± va/cs 



(52) 



where we have used rj = c^/47rcr||. If ions and electrons are 
the only charged species, then 



X = 



VA 



Re^ 



(l + /3|)V2(l + /3| 



2\l/2 



(53) 



In the resistive regime, |/3| <C 1 for both charged species 
and (53) shows that the criterion for non-Hall MRI pertur- 
bations to grow (W99), X > va/cs, is equivalent to Rcm > 1. 

We can now obtain an expression for x in terms of the 
electron fraction Xe = rie/nn at the midplane. 



X = 



1 B'ai. 



Pi-f3e 



eUeB 

(l+/3|)l/2(l+/32)l/2 

Ue < av > me 



(54) 



where we used expression (20) for cr± and assumed Pi -ti |/3e 
In the above equation 



< a-t; >w 1 X 10 



i28fer\^''^ 



(55) 



is the momentum-transfer rate coefficient for electron- 
neutral scattering. Note that the dependence of x on mag- 
netic field strength now enters only through the electron Hall 
parameter /3e. Following Fromang, Terquem & Balbus (2003) 
we assume that grains have settled out and the electron (and 
ion) number density is determined by the recombination of 
metal ions given by, 

1/2 

« I I , (56) 



C,nu 
a 



where a « 3 x 10 ^^^cm^s ^ is the radiative 
recombination rate for metal ions. The ionisation rate 
^ is assumed to be due to cosmic rays at a rate 
10~i^exp(— E/96gcm~2)s~i H~i where E is the disc sur- 
face density. This dominates x-ray ionisation for the column 
densities we shall consider here. Results are shown in Table 
1 for a nominal 1 solar mass star and B = 10 mG. 

For this strength of the magnetic field, |/3e| <g; 1 at 1 
AU and x scales as B^ (up to about 200 mG) at this radius. 
On the other hand, /3e is greater than 1 at 5 and 10 AU, so 
X will scale linearly with B (sec equation 54). In particular 
for B = 100 mG at 1 AU, X ~ 0.00088, consistent with the 
detailed calculations in Wardle (2003). 



4 BOUNDARY CONDITIONS 

To solve equations (39) to (42) it is necessary first to inte- 
grate the system of ordinary differential equations (ODE) in 
(39). This problem can be treated as a two-point boundary 
value problem for coupled ODE. Five boundary conditions 
must be formulated, prescribed either at the midplane or at 
the surface of the disc. 

At the midplane: A set of boundary conditions can 
be arrived at by asuming fluid variables have either 'odd' 
or 'even' symmetry about the midplane. 'Odd' symmetry 
means the variable is an odd function of z and vanishes at 
z = 0. Conversely, when 'even' symmetry is applied, the vari- 
able is assumed to be an even function of z and its gradient 
is zero at the midplane. In this paper we applied the odd - 
even symmetry criteria to the perturbations in the magnetic 
field, SB{z) = ±SB{~z), where the upper (lower) sign cor- 
responds to even (odd) symmetry conditions. This contrasts 
with Lovelace, Wang & Sulkancn (1987), who applied the 
symmetry criteria to the fiux function ^'(r, z) = rA^f, with 
Ati, the toroidal component of the vector potential and ob- 
tained dBr.d>i'i", z) ~ -fSBr.ip{r, — z) as their symmetry cri- 
teria. The symmetry of a particular fiuid variable is assigned 
arbitrarily, subject to the constraint that fluid equations are 
satisfied. This means that two sets of boundary conditions 
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Table 1. Comparison of the magnetic coupling parameter Xo and the ionisation fraction Xe at the midplane for different radial positions 
ro, with B = 10 mG and assuming grains have settled out (Promang, Terquem & Balbus 2002). Also shown axe the assumed temperature 
To and calculated values of nn, va/cs, Ci Pe and Pi. 



ro (AU) 


To (K) 


uh (cm 3) 




C (s-iH- 


') 




ft 




Xo 


1 


280 


6 X 10" 


0.00076 


5.76 X 10 


-22 


0.035 


7.68 X 10-"' 


7.32 X 10^13 


8.9 X 10-s 


5 


130 


7 X 10^2 


0.010 


4.81 X 10^ 


-18 


4.44 


6.58 X 10"^ 


5.11 X 10-1° 


1.9 


10 


90 


1 X 10i2 


0.033 


9.76 X 10 


-18 


37.31 


0.046 


1.76 X 10"^ 


19 



are equally valid, obtained by reversing the assumed sym- 
metry of the fluid variables. Perturbations obtained with a 
particular set of boundary conditions are displaced a quarter 
of a wavelength from those found with the other one. The 
growth rates of these solutions lie at intermediate points of 
the curve i' vs k obtained from the local analysis (W99), 
as expected. Evidently, no generality is lost by focusing in 
one of these two possible sets of solutions. We compare the 
growth rate versus number of nodes of perturbations ob- 
tained with 'odd' and 'even' symmetry in section 5.2 (Com- 
parison with local analysis). For the rest of the analysis pre- 
sented in this paper we chose to assign odd symmetry to 
SBr and 5B^, so they vanish at « = 0. This gives us two 
boundary conditions at the midplane. As the equations are 
linear, their overall scaling is arbitrary, so a third boundary 
condition can be obtained by setting one of the fluid vari- 
ables to any convenient value. To that effect, we assigned a 
value of 1 to SE'^.. Summarising, three boundary conditions 
are applied at the midplane: 

SBr = SBrj, = 0, and 

5K = 1. 

At the surface: At sufficiently high z above the mid- 
plane, ideal MHD conditions hold. This assumption is appro- 
priate in this case because the local coupling parameter x is 
inversely proportional to the density and so it is stronger at 
higher z regions where the density is smaller. When x > 10 
the growth rate and characteristic wavenumber of unstable 
modes differ little from the ideal limit (W99), so even though 
for simplicity we have assumed the conductivity tensor to be 
spatially constant, we can assume that flux- freezing condi- 
tions hold at the surface and the local dispersion relation 
is kvA = n (Balbus & Hawley 1991). As the Alfven speed 
increases with p'^^^, the wavelengths of magnetic field per- 
turbations increase with z and given the dependency of p 
with z (see equation 28), must tend to infinity as z — » oo . 
The displacements in the plane of the disc of an infinitely 
stretched perturbation should effectively vanish, so 5Br and 
5B^ should be zero at infinity. This gives us the remaining 
two boundary conditions required to integrate the system of 
equations (39). Consistently with them, both 5E' and 5 J 
vanish as well. 

Interestingly, this solution is consistent with 5E and Sv be- 
ing non-zero at infinity. The only requirement is that the 
gradient of the velocity in the vertical direction d5v/dz be 
zero when a — > oo, to prevent any horizontal stretching of 
the magnetic field. It may seem puzzling at first that 5v 
is non-vanishing at infinity. This can be understood taking 
into account that these perturbations travel to infinity in a 
finite time too given by. 



Table 2. Comparison of maximum growth rate Vmax and imm- 
ber of nodes A'^ of the fastest growing modes for all conductivity 
regimes and two different locations of the boundary. In all cases 
va/cs = 0.1 and Xo = 10. 





z/H = 5 




z/H = 7 




Conductivity Regime 


^max 


Af 




N 


Ambipolar Diffusion 


0.7303865 


5 


0.7303869 


5 


Hall limit (aiB^ > 0) 


0.7498761 


5 


0.7498761 


5 


Hall limit (aiB^ < 0) 


0.7461857 


5 


0.7461854 


5 


Comparable conductivities 


0.7345455 


5 


0.7345459 


5 


Opposite conductivities 


0.7340035 


5 


0.7340039 


5 



f°° dz H ( z'\ H 

too = / = — / exp — - \dz = Vtt — , 

^0 '<^Al VA Jo V 4 y VA 



(57) 



where z is the vertical coordinate in units of the scaleheight 
H and vai is the local value of the Alfven speed. Because of 
this finite travel time to infinity, the fiuid can retain a finite 
velocity when a — > oo. Furthermore, through equations (37) 
and (38) it is clear that 5E, the perturbations in the electric 
field as seen from the laboratory frame, are finite at infinity 
as well. 

These boundary conditions are strictly valid at infinity, but 
will also hold at a boundary located sufficiently high above 
the midplane. We chose to locate the boundary at z/H = 5 
after confirming that increasing this height does not sig- 
nificantly affect either the structure or the growth rate of 
unstable modes. This can be appreciated in Table 2, which 
compares the maximum growth rate Vmax and number of 
nodes N (proxy for wavenumber) of the perturbations in 
all conductivity regimes for two different locations of the 
boundary. Summarising, the boundary conditions adopted 
at the surface are: 

5Br = 5B4, = Q, 

at z/H = 5. 

This system of equations is solved as a two-point bound- 
ary value problem for coupled ODE by 'shooting' from the 
midplane to the surface of the disc and simultaneously ad- 
justing the growth rate v and 5E'^ until the solution con- 
verges. 



5 RESULTS 
5.1 Test Models 

We solved the system of equations (39) for different con- 
ductivity regimes, coupling between fluid components and 
initial magnetic field strengths. As discussed in section 3.1, 
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Table 3. Relative values of the components of the conductivity 
tensor perpendicular to the magnetic field ui and (T2 and fiducial 
values of the coupling parameters Xo (good coupling, poor cou- 
pling limits) and ua/cs adopted to explore the structure of the 
perturbations for all conductivity regimes. 



Conductivity Regime 


(T 






Xo 


va/cs 


Ambipolar Diffusion 


CTl 


= 




10, 1 


0.1 


Hall Limit aiB^ > 


0-2 


= 0, ai 


> 


10, 1 


0.1 


Hall Limit cri < 


CT2 


= 0, CTl 


< 


10, 2 


0.1 


Comparable conductivities 


(Tl 


= T2 




10, 1 


0.1 


Opposite conductivities 


(Tl 


= —ai 




10, 1 


0.1 



under the linear approximation and disc model adopted in 
the present paper, the ambipolar diffusion and resistive con- 
ductivity regimes are identical, so even though throughout 
this work wc have labelled the case when (Ti = as the 
'ambipolar diffusion' limit, it should be kept in mind that 
this condition describes the resistive regime as well. Two 
different Hall limits exist, as the growth rate of the MRI 
depends on the orientation of the initial magnetic field with 
respect to the disc angular velocity vector fl (W99). The 
case when Bo is parallel (antiparallcl) to fl is characterised 
by aiB, > {aiB^ < 0). 

We calculated the growth rate and vertical structure of 
all unstable perturbations for different conductivity regimes 
with va/cs = 0.1. The degree of coupling between the mag- 
netic field and the neutral component of the fluid was char- 
acterised by either Xo = 10 (good coupling) or Xo = 2 or 1 
(poor coupling). The choice of Xo for the low coupling anal- 
ysis is dependent on the conductivity regime of the fluid. 
We took Xo = 1 for all regimes, except the Hall ((J\Bz < 0) 
limit, where Xo = 2 was adopted as our code fails to con- 
verge for X < 2. In this regime, the local analysis shows 
that when 0.5 < x < 2 all wavenumbers grow (W99). We 
believe that this complex structure of the perturbations at 
ever increasing k prevents our code from converging when 
Xo<2. 

We also examined the dependency of the structure 

of the fastest growing modes, their growth rate and the 
height of maximum amplitude, with the coupling Xo and the 
strength of the magnetic field for all conductivity regimes. 
To study the effect of the magnetic coupling, the value of 
va/cs was fixed at 0.1. The impact of the strength of the 
field was explored for good and poor coupling conditions. 

Finally, we studied the dependency of the growth rate 
of the most unstable perturbations with the coupling Xo for 
va/cs = 0.1 and 0.01 and with the magnetic field strength 
for Xo = 10, 2 and 0.1. 

The relative values of aj and a2 used to characterise 
each conductivity regime, together with the values of va/cs 
and Xo used to explore the structure of the perturbations 
are summarised in Table 3. 
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Figure 1. Growth rate versus number of nodes (proxy for 
wavenumber) of the MRI for different conductivity regimes and 
coupling at the midplane Xo- Circles show the ambipolar diffusion 
limit {(71 = 0) and triangles the Hall limit ((T2 = 0, aiBz > 0). 
Filled symbols correspond to the good coupling case Xo = 10 and 
open ones to the poor coupling case Xo = 1- Crosses show the am- 
bipolar diffusion limit with 'even' boundary conditions applied to 
SBr and SB^ . Note that results in this ease are displaced a quar- 
ter of a wavelength (one node) from those obtained with 'odd' 
boundary conditions. In all cases va/cs = 0.1. 



to z — +5, as a proxy for the wavenumber k to compare 
our results with those of W99. Results are shown in Fig. 1 
for the ambipolar diffusion and Hall (aiB^ > 0) limits for 
good and poor coupling. In both regimes the reduction of 
the wavenumber of the fastest growing perturbation with Xo 
is obtained, as expected from the local analysis. Reducing 
Xo also diminishes the growth rate of the instability in the 
ambipolar diflfusion limit, but v remains unchanged for the 
Hall regime, as expected from the local results. 

These results confirm our expectation that applying 
boundary conditions and integrating the fluid equations in 
the vertical direction would restrict the unstable frequencies 
from the continuous curve v vs k obtained in the local anal- 
ysis to a discrete subset of global unstable modes supported 
by the fluid. 

Also shown as crosses in Fig. 1 are the growth rates 
obtained with 'even' boundary conditions applied to SBr 
and SB^ at the midplane (see section 4), for the ambipolar 
diffusion limit with Xo = 10. As expected, the perturbations 
are displaced a quarter of a wavelength (one node) from 
those obtained with 'odd' boundary conditions. 



5.2 Comparison with local analysis 

The linear growth of the MRI as a function of wavenumber 
in a local analysis shows that p vsk generally takes the form 
of inverted quadratics (W99). The local wavenumber of the 
perturbations change with z, so we use the number of nodes 
of 5Br over the entire thickness of the disc, from z = —h 



5.3 Structure of the Perturbations 

Fig. 2 shows the perturbations in all fluid variables as a func- 
tion of height, from the midplane [z/H = 0) to the surface 
of the disc (z/H = 5), for the ambipolar diffusion regime 
and good coupling (xo = 10). They are obtained through 
equations (40) - (42) once the ODE system (39) has been 
integrated. Note the non-zero values of 5v and 5E at the 
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Ambipolar Diffusion 



Hall CT.B >0 



Figure 2. Structure of the perturbations in all fluid variables as 
a function of height for the most unstable mode in the ambipolar 
diffusion limit, for good coupling (xo = 10) and va/cs = 0.1. 
The growth rate is u = 0.7304. Note the non-zero values of dv 
and SE' at the surface, due to the finite travel time to infinity of 
the perturbations. 



surface, as discussed in section 4. From this point onwards, 
discussion will be focussed in the perturbations of the mag- 
netic field only. Unless otherwise stated, the radial compo- 
nent of the field 5Br is plotted with a solid line and the 
azirnuthal component SB^ with a dashed lino. As the over- 
all scale of our linear equations is arbitrary, plots depicting 
the structure of the perturbations either do not show the 
scale of the vertical axis (corresponding to the amplitude of 
the perturbations) or show a conveniently normalised scale, 
for reference purposes. 



5.3.1 Effect of the Conductivity Regime 

Fig. 3 compares the structure of all unstable perturbations 
for the ambipolar diffusion and Hall {aiB^ > 0) regimes 
under good coupling {xo = 10). The modes are ordered by 
the number of nodes. We find that at small v, the struc- 
ture of the perturbations in both cases is very similar, but 
significant differences arise when the growth rate is close to 
maximum. Then, ambipolar diffusion perturbations peak at 
the node closest to the surface, while Hall {a^B^ > 0) ones 
peak closest to the midplane. This behaviour is linked to 
the change in the local coupling x with z and its effect in 
the structure of the perturbations for different conductivity 
regimes. It will be discussed further in the next section. At 
this good coupling level, there are no appreciable differences 
in the structure or growth rate of perturbations between 
both Hall limits results, as expected from the local analy- 
sis (W99) . When both conductivity components arc present 
(cri = ±(72 cases), the structure of the perturbations is sim- 
ilar to the ambipolar diffusion limit. This property is also 
dependent on the value of the local coupling x and will be 




, / 


7 \ 


u = 0.5153 

1 , , , , 1 < , , , 1 , < , < 1 < , < , 1 , 


/ 

A/ 
^ \ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

\ 


^- u = 0.6326 

1 . . . . 1 , . . . 1 . , . , 1 , . , . 1 . 


A 
,1,1, 


11 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
\. \ A 


\ v.' u = 0.6983 

iV , , , 1 1 , , , 1 , 1 , 1 1 1 , 1 , 1 , 


1 1 1 1 1 


1 1^1 1 Ml 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
A ■ \ / \ 





XI- 

\ \ ■ 1^ = 1). rou4 

1 , V, , i-/, 


1 1 1 1 1 
A f 

.i 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


— ^ 

V = 0.6908 

1 , , , , 1 < , , , 1 , < , < 1 < , < , 1 , 


n 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


, 1 ,v. 


y " V V = 0.6310 




HI 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

i iV\^ ^ 


1 V 1/ = 0.5098 






, 1 1 1 1 , 1 1 1 1 1 1 1 1 1 , 1 1 1 1 1 1 


1 , , , , 1 < , , , 1 , < , < 1 < , < , 1 , 


1 2 3 4 5 
z/H 





Figure 3. Structure and growth rate of all unstable modes of 
the MRI for the ambipolaj: diffusion (<ti = 0) and Hall (1T2 = 0, 
a^Bz > 0) cases. In all plots Xo = 10 and va/cs = 0.1. For 
this good coupling, there are no differences between both Hall 
limits. When both conductivity components are present {cri = 
±(T2 cases), results resemble the ambipolar diffusion limit shown. 
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analysed in the next section. We found 9 to 10 unstable 
perturbations in all cases. 

Under low coupling conditions, Xo = 2 or 1 depending 
on the conductivity regime (see Table 3), fewer unstable 
modes grow for both the ambipolar diffusion, Hall {ajBz > 
0) and the comparable conductivities (ui = a2) regimes. 
Six to eight unstable perturbations are found in these cases. 
As expected from the local analysis (W99), the range of 
wavenumbers for which unstable modes exist is reduced eis 
compared with the good coupling cases. 

Results are quite different for the two remaining con- 
ductivity regimes (Hall aiB^ < and ui = —172). More 
unstable modes arc found in these cases; 12 in the opposite 
conductivities case (cri = — (T2) and a total of 27 for the 
Hall {aiBz < 0) limit. In this last case in particular, despite 
the low coupling of ionised and neutral components of the 
fluid, unstable modes have a very complex structure (high 
wavenumber). Fig. 4 shows the structure of two such modes 
at the low growth rate, high wavcrmmbcr region of the v 
- k space. Note that unstable modes are so closely spaced 
that increasing the number of nodes by a few only magi- 
nally changes their wavenumber and has little effect on their 
growth rate. This complexity is expected from the form of 
the dispersion relation at low coupling from the local anal- 
ysis (W99). Non-linear simulations (Sano & Stone 2002a, 
b) confirm that the many growing modes in this regime 
strongly interact with each other and the instability devel- 
ops into MHD turbulence. This turbulence is a transient 
phase that eventually dies away in two-dimensional simula- 
tions (Sano & Stone 2002a) but it is sustained in full three- 
dimensional models (Sano & Stone 2002b). In both cases the 
non-emergence of the typical two-channel flow obtained in 
other regimes is noted by the authors. 

Finally, some of the perturbations show a structure re- 
sembling an interference pattern (see Fig. 5). They were ob- 
tained specifically in regimes where ambipolar diffusion is 
present, for good and poor coupling, but not in any of the 
Hall limits. This pattern can be explained recalling that lo- 
cal results show that two unstable modes exist with the same 
growth rate and different wavenumber. Despite this, just one 
global mode is found for each f in this analysis. Again, the 
application of boundary conditions and integration along 
the vertical direction restricts global unstable modes from 
those possible under a local analysis. The interference pat- 
tern suggests that global modes are a superposition of two 
WKB modes with (nearly) the same growth rate and which 
are not global solutions themselves. The interference of Fig. 
5 was successfully replicated through the superposition of 
two local modes with i> = 0.7004 using the analytical expre- 
sions in W99 for the ambipolar diffusion limit. 

5.3.2 Effect of the coupling parameter Xo 

Fig. 6 compares the structure and growth rate of the most 
unstable modes of the MRI for all conductivity regimes 
as a function of the coupling parameter Xo- In all cases 
va/cs = 0.1. We notice that reducing the coupling Xo 
causes the wavenumber (i.e. the number of nodes) of un- 
stable modes to diminish in all conductivity regimes except 
the Hall aiBz < limit (rightmost column of Fig. 6), for 
which this dependency is inverted. The growth rate is also 
reduced at a rate that depends on the conductivity regime 




Figure 4. Structure of two unstable modes in the Hall (uiBz < 
0) limit for poor coupling (xo = 2) and va/cs = 0.1. Note the 
complex structure of the perturbations (high wavenumber). Un- 
stable modes are so closely spaced that increasing the number 
of nodes only marginally changes their wavenumber and growth 
rate. 



of the fluid. These results are expected from the findings of 
the local analysis (W99). 

It is evident from Fig. 6 that at very high magnetic 
coupling (xo ~ 100), the fluid is close to ideal MHD con- 
ditions and results obtained in all conductivity regimes are 
alike. When the coupling is reduced to Xo ~ 10, we be- 
gin to appreciate differences between them. In particular, 
the amplitude of the perturbations when ambipolar diffu- 
sion is present peaks close to the surface while in both Hall 
limits the maximum amplitude is closest to the midplane. 
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Figure 5. Structure of an unstable mode in the ambipolar diffu- 
sion limit ((Ti = 0) for good coupling (xo = 10) and v^/cs = 0.01. 
Note the interference pattern of tlie perturbation, which suggests 
that this global mode is a superposition of two WKB modes of 
similar growth rates. 



This is more clearly appreciated in Fig. 7, which plots the 
height of maximum amplitude of the fastest growing modes 
as a function of the coupling parameter Xo and the con- 
ductivity regime of the fluid for va/cs = 0.1. This figure 
shows that pure Hall regimes (aiBz > and aiB^ < 0) 
peak closer to the midplanc, for all Xo studied, than the 
cases when ambipolar diffusion is present. This behaviour 
can be explained by the dependency of the local growth of 
the instability with x for different conductivity regimes. The 
maximum growth rate of ambipolar diffusion perturbations 
increases with the local x (W99), which in turn is a func- 
tion of height. As a result, at higher z, the local growth 
of the instability increases, driving the amplitude of global 
perturbations to increase. Hall {aiBz > 0) perturbations, 
on the contrary, have the same Vmax for all X: so the in- 
stability is not driven from any particular vertical location, 
which explains the flatness of their envelope. This also ex- 
plains why in Fig. 3 the differences between these regimes 
arc apparent at close to maximum growth: The increment 
in the local growth rate with the coupling x is less marked 
for slow growirrg perturbations. This is also appreciated in 
the local analysis by the form of the dispersion relation for 
different x in the ambipolar diffusion regime (W99). 

It is also clear from Fig. 7 that perturbations in the am- 
bipolar diffusion and Hall {aiBz > 0) limits peak at higher 
z when Xo is reduced. This dependency is driven by the re- 
duction of the wavenumber of most unstable perturbations 
with the coupling in these regimes (W99, see also Fig. 6). On 
the contrary, in the Hall (aiBz < 0) limit, the wavenumber 
of the fastest growing mode increases as Xo is reduced and 
the perturbations peak closer to the midplane with weaker 

Xo- 
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Figure 7. Height of maximum amplitude of the most unstable 
perturbations as a function of Xo for all conductivity regimes. In 
all cases va/cs = 0.1 



Finally, looking at the first throe columns of Fig. 6 
it is evident that the structure of unstable modes in the 
(Ti = ±(T2 conductivity regimes are remarkably similar to 
the ambipolar diffusion limit shown in the leftmost column 
of the figure for Xo ^ va/cs = 0.1. When the coupling is 
weaker than this value, the structure of unstable modes in 
these regimes is no longer alike. To explain this we recall 
that the minimum degree of coupling for unstable modes to 
grow, determined by the requirement that a wavelength fit 
in the disc scalelicight, is given by x ^ va/cs in the ambipo- 
lar diffusion limit, and x ^ va^ /cs^ in the Hall case (W99). 
For Xo < 0.1 then, the growth rate of ambipolar diffusion 
((Ti = 0) perturbations is expected to drop markedly and 
the envelope of the perturbations will bo mainly determined 
by the Hall effect. This transition in both ai = ±a2 cases, 
is clearly seen in Fig. 7. In both cases, for Xo ^ 0.1 the per- 
turbations resemble the ambipolar diffusion limit (<7i = 0). 
In the (71 = (72 regime, when Xo < va/cs the perturbations 
resemble those in the Hall {aiB^ > 0) limit (compare lowest 
panels in the second and fourth columns of Fig. 6), consis- 
tent with the notion that ambipolar diffusion effects are no 
longer important in this region of parameter space. This im- 
plies a constant i^max for weaker Xo, so the perturbations 
tend to peak closer to the midplane. On the contrary, in 
the CTi = —a2 case, Umax continues to diminish with x once 
X < va/cs (W99). In this case, the range of growing modes 
is always finite (as opposed to the Hall aiBz < limit) and 
there is a fastest growing mode for every x- As a result, the 
instability is driven at intermediate z (as in the ambipolar 
diffusion limit) and accordingly, the perturbations peat at 
higher z with decreasing Xo- The height of maximum ampli- 
tude increases faster as Xo is reduced in this regime than in 
the ambipolar diffusion limit. This occurs because Vmax for 
a given local x is greater in this case (W99, see also section 
5.4.1) so global perturbations are amplified even further. 
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Figure 6. Structure and growth rate of the most unstable modes of the MRI for all conductivity regimes and different values of Xo- In 
all cases va/cs = 0.1. The value of the coupling parameter Xo is indicated at the top right corner of each panel. The growth rate iy) of 
the perturbations is shown in the lower right corner. Note that the Hall aiBz < regime is explored for Xo > 2. In this conductivity 
regime our code fails to converge for x < 2 as in this region of parameter space the range of wavenumbers for which local unstable modes 
exist becomes infinite (W99). 
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5.3.3 Effect of the Magnetic Field Strength 

In ideal MHD, the weaker the magnetic field the higher the 
minimum wavenumber of the perturbations (Balbus & Haw- 
ley 1991). The results of this study are consistent with this 
finding. With va/cs ~ 0.005 the perturbations grow with 
very high wavenumbers in all conductivity regimes. This can 
be appreciated in Fig. 8 and 9 for the ambipolar diffusion 
and Hall {aiBz > 0) limits, respectively, for good coupling 
(Xo = 10). At this Xo, solutions for both Hall limits are sim- 
ilar. Also, (Ti = ±iT2 regimes are similar to the ambipolar 
diffusion limit, as expected. Note the interference pattern of 
perturbations in lower panels of Fig. 8. 

We also studied the dependency of the height of maxi- 
mum amplitude with the strength of the magnetic field (Fig. 
10) for good (xo = 10) and poor coupling (xo = 2). In the 
former case, when ambipolar diffusion is present, with and 
without the Hall effect and regardless of the sign of aiB^, 
the location of maximum amplitude of the fastest growing 
perturbations as a function of va/cs is similar, which is ex- 
pected as Xo > va/cs (see Fig. 10, panels a, b and c). We 
obtained unstable modes for va/cs up to 1. As va/cs is re- 
duced from this value, the perturbations peak at higher z 
until va/cs ~ 0.04. However, the location of maximum am- 
plitude begins to diminish as the field is further reduced, 
which could be caused by the interference pattern, and very 
high wavenumber, of the perturbations (see also Fig. 8). The 
height of maximum amplitude peaks at z ^ 2.8. 

In both Hall limits the height of maximum ampli- 
tude increases with the strength of the magnetic field until 
va/cs ~ 0.3 and then remains unaffected with further in- 
crements of va/cs (panels d and e of Fig. 10). This occurs 
because as the strength of the magnetic field increases, the 
wavenumber of the perturbations diminish, which pushes 
the maximum amplitude to higher z. For va/cs ^ 0.3 the 
perturbations have only one node, so any further increase in 
the magnetic field strength have little effect in the location 
of the max;imum amplitude. 

In the low coupling case results are very similar to the 
Xo ~ 10 cases. We note that in regimes where ambipolar 
diffusion is present (left hand side panels of Fig. 10), per- 
turbations tend to peak at a lower z/H than in the good 
coupling cases when va/cs < 0.1. In the Hall {aiBz < 0) 
regime, there are unstable modes for 'suprathermal' field 
strengths {va/cs up to 2.9). This will be further analysed 
in section 5.4.2, dealing with the dependency of the growth 
rate of this instability with the strength of the magnetic 
field. In this case, the height of maximum amplitude peaks 
at va/cs « 0.5 and then gradually diminishes as the field is 
incremented beyond this value. 



5.4 The perturbations in parameter space 

5.4.1 Effect of the coupling 

Fig. 11 shows the growth rate of the most unstable modes 
as a function of the coupling Xo for all conductivity regimes 
with va/cs = 0.1 (top panel) and 0.01 (bottom panel). The 
Hall {aiBz < 0) limit could not be modelled for Xo < 2 
because in this region of parameter space the range of 
wavenumbers for which unstable modes grow becomes in- 
finite (see section 5.1). We find that at the good coupling 
limit the instability grows at a rate similar to its ideal value 
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of 0.75Q for all conductivity regimes. As the coupling dimin- 
ishes, the growth rate is reduced at a rate that depends on 
the conductivity regime of the fluid. In the Hall {aiB^ > 0) 
case the growth rate remains unaffected until Xo ~ 0.01, and 
then diminishes drastically to O.IQ when Xo ~ 0.005. The 
ambipolar diffusion limit has a much more gradual reduction 
of h'max witli Xo- In this case, the growtli rate departs sig- 
nificantly from the ideal value for Xo ~ 0.1 and then drops 
rapidly, reaching ~ 0.007f2 for Xo ~ 0.008. This is in agree- 
ment with findings by W99 that unstable modes grow when 
X ^ va/cs in the ambipolar diffusion limit, and x ^ va' jcs^ 
in the Hall case (see also section 5.3.2). When x is less than 
these values, perturbations are strongly damped. 

When va/cs ~ 0.01 (bottom panel of Fig. 11), the 
growth of Hall (aiBz > 0) perturbations is constant at 
about the ideal rate 0.75n until Xo ~ lO"**. Below this, it 
plummets to zero as expected. On the contrary, in the am- 
bipolar diffusion and a\ — ±(72 regimes, Vmax begins to di- 
minish much sooner. It is also noticed that, in the ambipolar 
diffusion and ai = ai cases, the growth rate increases again 
after reaching a minimum for Xo ~ 0.05. This is caused 
by the high wavenumber of the perturbations due to the 
weakness of the magnetic field. In these conditions global 
effects (stratification) arc less important and the maximum 
growth rate diminishes with the coupling as per local results 
(W99). As the wavenumbers of the perturbations decrease 
when Xo is reduced, for sufficiently low Xo ( ~ 0.05), fe is low 
enough for global effects to be important again and modify 
the growth rate of unstable perturbations. With the local 
X increasing with height above the midplane, stratification 
will tend to increase the growth of global modes at low Xo- 



-0.5 



o 



-1.5 



A 


1 1 1 1 1 1 1 1 1 1 1 1 
A A A A A ^ _0 S § # " 
S * 

□ ° 8 


A 

7 A ° 


a ^ V /c =0 1 " 
* A/ s ■ 








A 
13 


*c7j = L 


□ 




□ 


A cr2=0, (7iB^>0 \ 


* 




□ 


oCT2 = 0, (7iB^<0 7 


* 








* 

1 1 


1 1 1 1 1 1 1 1 1 1 1 1 


-2 


-1 1 



-0.2 
-0.4 



^ -0.6 

O 



-0.8 



-1.2 



A A 



A^AA AAA AAA AA O© © 



* 



Va/c,=0.01 



-4 



-2 

Log(Xo) 



Figure 11. Growtli rate of the fastest growing modes of the MRI 
as a function of Xo for different conductivity regimes, va/cs =0.1 
(top panel) and 0.01 (bottom panel). 



5.4.-2 Effect of the magnetic field strength 

The dependency of the maximum growth rate with the 

strength of the magnetic field for all conductivity regimes 
is shown in Fig. 12 for Xo = 10, 2 and 0.1. In the good cou- 
pling case (top panel), increasing the strength of the mag- 
netic field has little effect in the growth rate of the most 
unstable modes of all conductivity regimes until va/cs ~ 1, 
where it drastically drops to zero. These results are similar 
to the ideal MHD case, which predicts that at this strength 
of the magnetic field the wavelength of most unstable modes 
become ~ H, the scaleheight of the disc, and the perturba- 
tions are strongly damped. 

In the Xo = 2 case shown in the middle panel, we found 
unstable modes in the Hall limit (aiBz < 0) for va/cs up to 
2.9. We know from the local analysis (W99) that once the 
local X ^ 2, unstable modes exist for every kvA/^i in this 
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conductivity regime. As a result, even for suprathermal fields 
{va/cs > 1), there are still unstable modes with kH < 1 
growing within the disc. 

Results for Xo = 0-'^ (bottom panel) show clearly how 
i^max plummets when va/cs > x (ambipolar diffusion limit) 
or va' jcs^ > X (Hall limit), £is expected. 

Finally, turning our attention to the dependency of the 
growth rate of the instability with the field strength at low 
va/cs, wc appreciate in all panels of Fig. 12 that Vmax ini- 
tially increases as 'ua/cs is reduced, until it reaches a maxi- 
mum. Further reductions in va /cs cause the growth rate to 
diminish monotonically. Comparing these results with the 
findings of the local analysis (W99), we can show that the 
growth rates of global unstable modes at weak magnetic 
field strengths tend to the local values for the relevant cou- 
pling Xo- This can be explained simply by the increase in 
wavenumber of unstable modes as va/cb is reduced, which 
causes global effects (driven by stratification) to be less im- 
portant. As a result, the growth of global modes does not 
differ significantly from the local values at the same cou- 
pUng. 



6 DISCUSSION 

The solutions presented in the previous sections illustrate 
the growth and structure of the MRI when diff'erent com- 
ponents of the conductivity tensor are dominant through- 
out the entire cross-section of the disc. Density stratification 
causes the local growth of unstable modes, and the ampli- 
tude of global perturbations, to be a function of height. The 
envelopes of short-wavelength solutions are shaped by this 
competition between different growth rates acting at dif- 
ferent vetical locations. Moreover, when x is weak at the 
midplane, long wavelength perturbations are important and 
vertical stratification is crucial in determining the growth of 
global MRI solutions. The results presented in this contri- 
bution confirm these expectations. 

When the Hall and Pedersen components are compara- 
ble, the Hall effect alters the structure and growth of unsta- 
ble perturbations for Xo < va/cs. In this region of parame- 
ter space, ambipolar diffusion perturbations have negligible 
growth, but unstable modes that include Hall conductivity 
still grow at v = 0.2 — 0.3. Hall aiB^ > perturbations 
grow faster than aiB^ < ones. Furthermore, under a weak 
magnetic field {va/cs < 0.01), the Hall effect significantly in- 
creases the growth rate of unstable modes at low coupling. 
When it dominates, unstable modes grow at close to the 
ideal rate for Xo ~ lO"''. 

The height above the midplane where the most unstable 
perturbations peak is dependent on the conductivity regime 
of the fluid. Consequently, the vertical location of the ac- 
tive zones within the disc, in which the MRI produces an- 
gular momentum transport and disc material is being ac- 
creted, is dependent on the configuration of the conductiv- 
ity tensor. Perturbations including ambipolar diffusion peak 
consistently higher than those in both Hall limits. Also, for 
Xo < va/cs, o\ = ±(72 modes peak at different heights, sig- 
nalling that the Hall cff'ect is dependent on the orientation of 
the magnetic field with respect to the disc angular velocity 
vector n. In this region of parameter space, when a^Bz < 0, 
modes peak at a higher z than when a^B^ > 0. The Hall 



effect does not significantly modify the dependency of the 
height of maximum amplitude of unstable modes with the 
strength of the magnetic field for good coupling. In this re- 
gion of parameter space ambipolar diffusion dominates and 
causes the perturbations to peak at higher z when va/cs is 
reduced. 

When aiBz < 0, Hall perturbations can have a very 
complex structure (high wavenumber), even at low coupling, 
and many modes are found to grow. In the non-linear stage, 
the interaction between these modes causes the MRI to de- 
velop into MHD turbulence with non-emergence of the typi- 
cal two-channel flow obtained in other conductivity regimes 
(Sano & Stone 2002a, 2002b). 

When does Hall diffusion determine the behaviour of 
the instability? Naively one might propose |cri|>(j2 as a 
criterion, but this ignores the level of coupling between the 
magnetic fleld and the gas. For example if x ^ 10, ideal MHD 
almost holds and there is little dependence of growth rate 
and structure on the diffusion regime (see Fig. 3). A use- 
ful criterion can be derived using the results of the local 
analysis in W99 and comparing the maximum growth rates 
with and without Hall diffusion in the weak-coupling limit. 
In the absence of Hall diffusion (i.e. ui =0), the maxi- 
mum growth rate for X ^ 1 is 3x/4. When Hall diffusion is 
present and |ci|/c±, the max:imum growth rate is either 
||ail/((Tx + 0-2) if criBz > or the instability is suppressed 
if aiBz < 0. In either case. Hall diffusion dominates the 
behaviour of the instability when 

X < Xcrit = ^ . (58) 

Thus even when IctiI <JC 0-2, the structure and growth rate 
of the magnetorotational instability are dominated by Hall 
diffusion if x < Ifil/crx. This is easily satisfied, for example, 
for the nominal conditions at the disc midplane 1 AU from 
the central protostar (sec Table 1), where |(Ti|/crx « |/3e| = 
0.035 (using eqs. 18 and 20 with |/3e| > 1) and x = 9 x 10"'^. 
At 5 AU, X ~ 2 and ai S> (72, so Hall diffusion dominates 
here also. Although these results depend on the assumed 
magnetic field strength, the conditions under which x < 
Xcrit arc so broad that wc can conclude that Hall diffusion 
determines the growth rate and structure of the instability 
over a large range of radii. 

Despite this. Hall diffusion has generally been neglected 
in studies of accretion discs in favour of the ambipolar diffu- 
sion or resistive limits. Here we illustrate the severity of this 
approximation by comparing the structure and growth rate 
of the unstable perturbations for a model with ai = 02 to 
those for 'simplified' pure ambipolar/resistive diffusion and 
Hall ((71 Bz > 0) models obtained by setting 02 or a\ to zero 
respectively and reducing the coupling parameter Xo by a 
factor of %/2 to reflect its dependence on (ai^ -f a-i'Y^^ (see 
eq. 50). The full model has va/c^ = 0.01 and Xo = 0.01414, 
whereas the the corresponding ambipolar diffusion and Hall 
limits have Xo = 'va/cs = 0.01 and the appropriate values 
of ai and (T2 as per Table 3. 

The comparison is presented in Fig. 13. The struc- 
ture and growth rate of pure ambipolar diffusion and Hall 
(fi-Bz > 0) perturbations are as expected for Xo ~ va/cs, 
with the ambipolar diffusion decaying towards the midplane 
and the envelope of the pure-Hall solutions being fairly con- 
stant. In both the ambipolar diffusion and the ai = a2 cases 
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we obtain a magnetic 'dead zone' near the midplanc where 
no perturbations grow (Gammie 1996, Wardle 1997). Com- 
paring the top panel of Fig. 13 with the middle and bottom 
ones it is clear that the Hall effect modifies both the struc- 
ture and growth of unstable modes. In particular, the extent 
of the dead zone is reduced and the growth rate is increased. 
According to these results, both the depth of the active zones 
within the disc and the rate of angular momentum transport 
by unstable modes can be significantly modified by the Hall 
effect. 

Although our solutions incorporate the effect of den- 
sity stratification on the coupling parameter and the Alfven 
speed, we have assumed that the components of the conduc- 
tivity tensor do not vary with height. While this simplifica- 
tion permitted us to compare the behaviour of the instability 
in different regimes, the conductivities in a real disc will re- 
fiect the height-dependence of charged particle abundances 
and their Hall parameters (see §2.1). Different regimes are 
expected to dominate at different heights (Wardle 2003). 

Nonetheless, it is clear from the simplified comparison 
presented here that Hall diffusion is an essential part of 
accretion in low-conductivity discs, and that it determines 
the extent of the magnetically-inactive 'dead zone' (Gam- 
mie 1996, Wardle 1997). Further, Hall diffusion will mod- 
ify any angular momentum transport within the dead zone 
that occurs via non-ajcisymmetric density waves driven by 
the active surface layers (Stone & Fleming 2003) because 
it will dominate the marginally magnetically-active regions 
of the disc just above the dead zone. Hall diffusion may 
therefore affect the ability of dust grains to settle towards 
the midplane and begin to assemble into planetesimals (e.g. 
Weidenschilhng & Cuzzi 1993). 



7 SUMMARY 

In this paper we have examined the structure and linear 
growth of the magnetorotational instability (MRI) in weakly 
ionised, stratified accretion discs, assuming an initially ver- 
tical magnetic field. This work is relevant for the study of 
low-conductivity accretion systems, such as protostellar and 
quiescent dwarf novae discs, where non-ideal MHD effects 
are important (Gammie & Menou 1998; Menou 2000; Stone 
et al. 2000). The formulation allows for a height-dependent 
conductivity but in this initial study we assumed the compo- 
nents of the conductivity tensor were constant with height. 
The analysis was restricted to perturbations with a vertical 
wavevector {k — kz), which arc the most unstable modes 
when initiated from a vertically aligned magnetic field (Bal- 
bus & Hawley 1991; Sano & Miyama 1999). In this case, the 
field-parallel component of the conductivity tensor plays no 
role and the ambipolar diffusion and resistive limits are iden- 
tical. The linearised system of ODE was integrated from the 
midplane to the surface of the disc under appropriate bound- 
ary conditions and global unstable modes were obtained. 
The parameters that control the evolution of the fluid are: 
{i) The coupling between ionised and neutral components of 
the fluid evaluated at the midplanc (xo), which relates the 
frequency at which non-ideal effects are important with the 
dynamical (keplerian) frequency of the disc; {ii) the mag- 
netic field strength characterised by the ratio va/cs at the 



midplane; and (in) the ratio of the components of the con- 
ductivity tensor perpendicular to the magnetic field a\/a2- 

In order to explore the growth and structure of unstable 
modes when different conductivity regimes dominate over 
the entire cross-section of the disc, we examined the follow- 
ing configurations of the conductivity tensor: ui = (the 
ambipolar diffusion or resistive limits), cr2 = (both Hall 
limits a^Bz > and aiB^ < 0), and the cases where both 
effects are important (cri = ±172). 

The main results of this study are highlighted below: 

(i) Global modes are a discrete subset from the contin- 
uous curve of possible vs. k combinations obtained with 
a local analysis (W99). These unstable modes can be ex- 
pressed as a superposition of two WKB modes with similar 
growth rate, which explains the interference patterns found 
in some of the perturbations. 

(ii) Ambipolar diffusion perturbations peak consistently 
higher above the midplane than solutions where Hall con- 
ductivity dominates. 

(iii) For good coupling (xo > va/cs), the structure and 
growth of the perturbations are mainly determined by am- 
bipolar diffusion. For a weaker coupling. Hall conductivity 
significantly modifies unstable modes. In this case, ai = (J2 
perturbations resemble the Hall aiBz > limit and peak 
closer to the midplane while ai = —a2 modes have their 
maximum amplitude closer to the surface. 

(iv) Hall limit [aiB^ < 0) perturbations can have a com- 
plex structure (high wavenumber) even for poor coupling 
(Xo = 2). There are also many unstable modes, which sup- 
ports findings that in this case the MRI evolves into MHD 
turbulence with non-emergence of the two-channel flow ob- 
tained in other regimes (Sano & Stone 2002a, 2002b). 

(v) As the coupling parameter Xo is reduced, departure 
from ideal growth 1/ ~ 0.75 occurs at a rate that depends 
on the conductivity regime. Hall limit perturbations grow 
at close to the ideal limit for Xo > va^/cs^- In the ambipo- 
lar diffusion approximation the growth rate decreases when 
Xo ^ va/cs. These results axe in agreement with predictions 
from W99. 

(vi) The weaker the magnetic field the higher the per- 
turbations peak in all regimes where ambipolar diffusion is 
present. On the contrary, both Hall limits peak closer to the 
surface with weaker va/Cs- 

(vii) Unstable modes grow when va/cs is increased until 
a critical value va/cs^^^^ is reached. At the critical va/cs 
the growth rate abruptly drops to zero. At good coupling 
VA/cs^^ff ~ 1 for all conductivity regimes. At the poor cou- 
pling limit (xo ~ 2), results are different only for the Hall 
regime {o\Bz < 0). In this case we obtain unstable modes 
for va/cs ~ 2.9. 

(viii) At very weak magnetic flelds {va/cs < 0.01), global 
effects are less important, due to the high wavenumber of the 
perturbations. In this region of parameter space the growth 
rates of MRI perturbations tend to the corresponding local 
values for the relevant fluid parameters. 

(ix) Hall difussion determines the growth of the MRI 
when x^ kil/'^'x- This condition is satisfied over a large 
range of radii in protostellar discs. 

(x) When the Hall regime dominates near the midplane 
and ambipolar diffusion is dominant closer to the surface, a 
larger section of the disc is unstable to MRI perturbations 
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and unstable modes grow faster than those obtained using 
the ambipolar diffusion approximation 
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Figure 12. Growth rate of the fastest growing modes of the MRI 
as a function of Vj^/cs for different conduetivity regimes. Prom top 
to bottom the coupling Xo is 10, 2 and 0.1. 
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Figure 13. Comparison of tlie structure and growth rate of the 
MRI for different configurations of the conductivity tensor for 
va/cs = 0.01. Top panel shows the case where both ambipolar 
diffusion and Hall {uiB^ > 0) conductivity terms are important. 
In this configuration the Hall regime is dominant close to the mid- 
plane and ambipolar diffusion dominates near the surface. Middle 
and bottom panels show the instability under the ambipolar dif- 
fusion and Hall approximations, respectively. In the top panel 



